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Abstract 

We construct analytic solutions of the low energy (i.e. tree-level) string cosmological effective 

action. We work with the "minimal" field content (i.e. graviton and dilaton) in the absence of 

^ ■ any dilaton potential and of any antisymmetric tensor field. Provided the metric is sufficiently in- 
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homogeneous we obtain solutions whose curvature invariants are bounded and everywhere defined 
in time and space. The dilaton coupling and its associated energy density are regular and homo- 

^ . geneous. A phase of growing (and non-singular) dilaton coupling is compatible with the regularity 

Q_i' of the curvature invariants without the addition of higher curvature (or higher genus) corrections 

^ . to the tree-level effective action. We discuss the symmetries and the physical properties of the 



X 



obtained solutions. 
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1 Introduction and motivations 

A quite intriguing aspect of cosmological models based on the low energy string theory effective 
action is the occurrence of a phase of growing dilaton coupling that might trigger a number of 
interesting implications which were specifically addressed in the framework of the pre-big-bang 
inflationary models |l|| . One of the distinctive aspects of this class of models is the amplification, 
via gravitational instability, of the vacuum inhomogeneities associated both with the metric and 
with the Abelian gauge fields [^. If this picture is correct, we should live, at later times in an 
Universe filled with stochastic backgrounds of gravitational waves and of of Abelian gauge fields 
which, thanks to the magnetic flux conservation in a hot plasma, might also have interesting 
effects. 

A lot of efforts have certainly been devoted to the analysis of the cosmological solutions of the 
low energy string theory effective action |^ in order to understand if the physical properties of 
those solutions were at all relevant for the explanation of some general features of our Universe. 
In all these studies one of the key assumptions was the complete homogeneity of the classical 
evolution of the geometry. It has been recently realized that it would be quite interesting to relax 
this very stringent assumption and, therefore, inhomogeneous string cosmological models were 
discussed from different points of view [^ |^ but always in the framework of the low energy string 
theory effective action. There are of course different motivations for these investigations. In gen- 
eral, it is always quite dangerous to make strict assumptions about the form of the cosmological 
solutions. We deal here with solutions of a completely non-linear system of evolution equations 
and we cannot exclude, in principle, physical situations were the Universe was very inhomoge- 
neous at early times evolving, ultimately, towards a homogeneous (and possibly isotropic) epoch 
at later times. In order to address these questions one should start, from the very beginning, with 
a completely inhomogeneous metric. Completely inhomogeneous configurations are also required 
in the analysis of the initial conditions of cosmological (inflationary) models. To relax the homo- 
geneity assumption turned out to be relevant for understanding if (and when) inflation could take 
place starting from more generic initial conditions. An example in this direction are the inves- 
tigations concerning the role played by initial (spherical) inhomogeneities in chaotic inflationary 
models |^. With the same motivations (but in the context of string inspired cosmological mod- 



els) quasi-homogenous metrics were used for the discussion of the initial value problem |^ in the 
pre-big-bang scenario |I|, 0]. In this last class of models, the discussion of the initial conditions in 
a completely homogeneous metric can lead to quite different results ^. We argue, therefore, that 
there is an urgent need of a better understanding of completely inhomogeneous configurations of 
the massless string modes. 

We cannot a priori exclude that the features of inhomogeneous string cosmological models 
might be significantly different from the ones of their completely homogeneous counterpart. 

An open question of (homogeneous) string cosmological models is the occurrence of curvature 
singularities [0 which turned out to be closely connected with the problem of a graceful exit 
from a phase of growing dilaton coupling 0, towards a phase of declining dilaton coupling. In 
a completely homogeneous metrics quite general (numerical) studies suggested that a phase of 
growing coupling is compatible with the regularity of the curvature invariants (in the String 
frame) provided higher order (curvature) corrections to the low energy (tree-level) effective action 
are included in the game 0. This idea received recently new attention ||10[. Moreover, the idea 
of a graceful exit transition from the initial phase of growing coupling to the subsequent standard 
radiation-dominated evolution has been the subject of many investigations [|TT] . 



In this paper we will not address the problem of the graceful exit but we will be concerned 
with the regularity of the curvature invariants. 

In the case of a completely inhomogeneous metric is it still mandatory to hit a singularity after 
a phase of growing coupling? Is it possible to regularize the curvature invariants if the Universe 
was sufficiently inhomogeneous at early times? 

These are the main question we will try to address and, in this sense, our approach is quali- 
tatively different from the ones followed in Ref. 0, |ll| (where homogeneity was assumed) and 
could be related to some pioneering attempts discussed in Ref. ||12|]. Regular (inhomogeneous) 
solutions of the low energy string theory effective action were constructed |12[ by "boosting away" 



the curvature singularities through the 0{d,d) symmetries which were extensively studied in the 
past [[131 , mi and which represent a natural generalization of the scale factor duality [jl], |16| . The 



key observation was that starting with a generic singular cosmology in 1 -|- 1 dimensions or with 
a 1 + 1 dimensional black hole, the singularity gets "boosted away" by 0(2,2) transformation 
involving a third (originally flat) dimension. Moreover, by considering the inhomogeneous cosmo- 



logical model of Nappi and Witten |jT5| it was also shown that the singularities of the model can 



be "boosted away" by a 0(3, 3) transformations involving a fifth dimension. The main conclusion 



of Ref. [0 was that, provided the coupling oscillates and provided an antisymmetric tensor back- 
ground is included, the curvature invariants are regular. Therefore, there are no doubts that it is 
certainly physically relevant to relax the homogeneity assumption usually made in string inspired 
cosmological models. Following this logical line and applying the Occam razor we want to study 
inhomogeneous string cosmological models with the minimal field content (i.e. the gravi-dilaton 
action). In other words, we want to understand if it is possible to have simultaneously satisfied 
these requirements: growing coupling, regularity of the curvature invariants, absence of antisym- 
metric tensor field, absence of any dilaton potential. In this sense the aim of our investigation is 



different from (but complementary to) the one of Ref. [|T^]. 

To tackle this problem many different strategies can be certainly employed. The most general 
one would consist in trying to solve the low energy beta functions for a generic metric g^i,{x,t) 
depending simultaneously upon space and time. A simplification adopted in previous studies [§, Q 
was also to stick to a diagonal inhomogeneous metric. Also the form of the inhomogeneous metric 
is, in a way, "minimal" . In fact the simplest non-trivial case of fully inhomogeneous metrics is 
the one where the inhomogeneities are distributed along a line. In this case the whole spatial 
dependence is reduced to only one coordinate (which we took, in our examples, to run along the 
the X axis). Needless to say that this is not the most general situation. At the same time there are 
no reasons to forbid this choice. A posteriori we finally found useful to impose an extra symmetry 
on the metric which will turn out to possess, in our discussion, two killing vector fields which 
are hypersurface orthogonal and orthogonal to each other. This form of the metric seems to be 
widely used in the study of inhomogeneous cosmological models in general relativity [|l7l. It is 
interesting to notice (as we will explicitly show) that this ansatz for the metric is compatible with 
a completely homogenous dilaton coupling which turns out to be, ultimately, only dependent upon 
time. 

We solved exactly the low energy beta functions in the String frame by using all the assumptions 
listed above. We computed the curvature invariants associated with the obtained solutions and 
we found a class of regular (and parity-invariant) solutions where the curvature invariants are 
all bounded. These solutions describe the physical situation where a growing dilaton coupling 



smoothly evolves between two asymptotically constant regimes. The physical parameter describing 
the solutions is essentially the maximal energy density of the dilaton background in string units. 

From a technical point of view, we have to mention that the most difficult part of our present 
exercise is to solve the off-diagonal component of the beta functions which do not contain second 
order derivatives with respect to time and which is then a constraint mixing first order time 
derivatives with first order spatial derivatives. Therefore, the first step will be the solution of 
the constraint. If the obtained solution will be compatible with the dilaton equation and with 
the remaining components of the beta functions, then, we will have a particular solution whose 
singularity properties can be investigated by the direct calculation of the curvature invariants. 

Before closing this introduction we want to recall that the work presented in this paper is cer- 
tainly overlapping with previous work performed in the context of general relativity. In particular 
in the limit of constant dilaton field our solutions match the well known inhomogeneous vacuum 
solutions of general relativity. The simplest inhomogeneous vacuum models are those with two 



space-like commuting killing vectors, known as orthogonally transitive G2 cosmologies ||T8|, [19 
Moreover, useful background material for the present investigation can be also found in [20| where 
soliton solutions were discussed in space-times with two space-like Killing fields. 

The plan of our paper is then the following. In Section 2 we will derive the explicit form of the 
evolution equations of the dilaton and of the metric in the String frame picture. In Section 3 we will 
briefly derive the dilaton vacuum solutions in the String frame and we will generalise them to the 
case of homogeneous dilaton background with inhomogeneous metric. In Section 4 we will focus 
our attention on the case of growing dilaton solutions and we will provide two physical examples 
whose main physical aspect is the parity invariance of the background geometry. Motivated the 
results of Section 4 we generalized our examples to a class of parity even solutions (Section 5) 
and we computed the curvature invariants in this situation. We found that they are regular as 
in the case of the examples of Section 4. Section 6 contains our concluding remarks. Because of 
the excessive length of the formulas contained in this paper we made the choice of reporting few 
(indeed relevant) technical results in the Appendix which then collects useful calculations for the 
interested reader. 



2 Basic Equations 



The low energy string theory effective action can be written in the String frame and in the absence 
of antisymmetric tensor as B: 
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R + g''^d^(f)dp(j) 



(2.1) 



where A^ is the string scale. We assume that the dilaton potential is completely negligible. We 
shall consider only the particular case of critical superstring theory with vanishing cosmological 
constant and six (frozen) internal dimensions. The evolution equations for the massless modes 
can be easily derived by varying the action with respect to the dilaton field and the metric: 



R - g'^^dacpdpcp + 200 = 0, 

Rf,u - ^QfivR + ^gi,vg'^^da(t)dii(t) + V^V,.^ - g^^v^ct) = 



(2.2) 
(2.3) 



(Vq, is the covariant derivative with respect to the metric (7^,^; Vq,V/30 = dadi3(f) — T'^pda(p; F^^ are 
the Christoffel symbols). Using Eq. ( p.2|) into Eq. ( p. 3D we find the familiar form of the tree- level 
beta functions 



R - g"'^d^(f)df3(f) + 200 = 0, 



r; + v,v'<p = o. 



(2.4) 
(2.5) 



We now specialize to the case of a diagonal metric which admits two commuting space-like killing 

B{x,t) 



vector fields, both of which are hypersurface orthogonal |T^ 

goo = A{x,t), g,cx = -A{x,t), gyy = -B{x,t)C{x,t), g. 
The corresponding line element is then 



C{x,t)' 



(2.6) 



ds^ = A{x,t)[dt^ - dx^] - B{x,t) 



C{x,t)dy^ + 



dz^ 



C{x,t) 



(2.7) 



From Eqs. ( p.6|) and (|2.7| ) the Christoffel symbols, the Riemann tensors and the Weyl tensors 
can be easily computed. In Appendix A we report the Christoffel symbols, the Ricci tensors 
and the curvature scalar which are required in order to express directly Eqs. ( p. 41 ) and ( p. 51 ) 
in the metric ( |2.6| ). In the metric ( p.6| ), Eqs ( p.4| ) and ( |2.5[ ) will produce a non-linear system of 



(partial) differential equations which we are going to solve. First of all suppose that the space-time 
dependence of the metric functions can be factorised: 

A{x,t) = a{t) a{x), B{x,t) = b{t) f3{x), C{x,t) = c{t) k{x), (2.8) 

(notice that the latin letters will denote the temporal part whereas the Greek ones will denote the 
spatial part of the metric functions). This decomposition has the great advantage, once inserted 
in Eqs. ( p.4| ) and ( |2.5| ), of transforming the original system into a system of non-linear ordinary 
differential equations in the two variables x and t. By inserting Eq ( p.8|) into the components of 
the Ricci tensors and into the curvature scalar given in Appendix A (see Eqs. ( |A.1| )-( |AI^ )) we can 
write down explicitly Eq. ( |2.4D and all the components of Eq. ( ^.5[ ) in the metric given in Eq. 
]6|). The (xx), iyy), (zz) and (00) components of Eq. (|275|) will be respectively: 
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whereas for Eq. <^A) and for the mixed component (i.e. (xO)) of Eq. (2.5) we 


will have 


respec- 


tively: 
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(in these equations and in all the paper we adopted the compact notation '■ ' = -§^ and ■ = ^)- In 
the system of Eqs. ( p.9D -( ^T3|) at the left hand side of each equation we only have tzme-dependent 



quantities whereas at the right hand side we only have s;?ace-dependent quantities. The exception 
is represented by Eq. ( |2.14| ) which is not a dynamical equation but a constraint mixing spatial 



and temporal derivatives. The coupled system of non-linear differential equations is written in 
the case where the dilaton coupling is completely homogeneous, namely 0' = 0, (t>{x-,t) = (f){t). 
The solution of the system of coupled, non-linear differential equations derived in the previous 
equations will be the main goal of the next Section. 

3 Homogeneous dilaton solutions 



The first step in the derivation of our solutions with homogeneous dilaton will be to consider 
the constant dilaton case where 0(t) = 0. In order to find a general expression able to solve the 
constraint, we sum and subtract Eqs. (p.lO| ) and ( p.llD , obtaining, respectively. 
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A particular solution of Eq. (|3.1| ) which might solve the constraint is 

b{t) = cosh /it, P{x) = sinhfix. 



(3.1) 



(3.2) 



It is true that Eq. (p.l|) can be solved also for another choices of trial functions: not all the 
possible choices are allowed by the off-diagonal constraint given in Eq. ( |2.14|) , which can now be 
re-written as 

2 tanh fit 
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Citanhyut, and 



fi 



tanh fix 



K 
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0. 



C2 tanh /it. 



(3.3) 



(3.4) 



the time dependence appearing in Eq. (3^) factorizes leading, ultimately, to a first-order (ordi- 
nary) differential relation involving only functions depending upon the spatial coordinates. There- 
fore direct (trivial) integration of Eqs. ( p.l|) and ( p.3|) , together with the conditions (|3.2| ) and (|3.4| ), 
gives us the following first integrals 



a 



C2 + C3 f. , fJ-X 

tanh — 
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H tanh — 

2 V 2 



ci + C2C3 - fici + fi^ 
2^ 



( tanh — 
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(3.5) 
(3.6) 



where C3 is now an integration constant and Ci and Ci are nothing but arbitrary constants which 
have to be fixed by requiring the compatibihty of Eqs. ( |3.2D and ( p.4| ) with the complete system 
given by by Eqs. ( p.9| )-( prT3D . Let us therefore find the compatibihty relations which have to be 



satisfied by Ci, c^ and C3. By direct substitution of Eqs. (|3.2|), (|3.4| ), ( p.5|) and (|3^ ) into Eqs. 
( ^I9D and (|2.12|) we get, respectively 

C3 = ?)^ + cl — 2/iCi, cl = 2/iCi + yU^. (3.7) 

Surprisingly enough the same two conditions also solve Eq. ( p.l3| ). Moreover Eqs. ( |2.10| ), ( p.ll| ) 
and ( p.l4| ) are automatically satisfied for every ci, C2 and C3 not leading (as expected) to any 
further relation among the constants appearing in our trial solution. The present solutions are 
the String frame version of analogous inhomogeneous vacuum solutions previously discussed in 
the context of general relativity |T^, |TB 



After having solved the case where the dilaton was constant we can now discuss the case where 
the dilaton is time dependent. The ansatz given in Eqs. ( p.2|) and ( p.4|) can be slightly modified 
in order to match the time dependence of the dilaton field, which now appears, in the equations 
of motion (|2.12|) - (|2.13|) with its first and second derivatives. The idea is then to assume that in 
the presence of the dilaton field the solutions of our system of equations can be written as 

f-j =citanh fit + T{t), [ - j = /itanh/xt + A(t), ( - J = catanh/xt + S(t) , (3.8) 

where, at the moment T{t), A(t) and S(t), are three unknown functions which should be hopefully 
determined by solving explicitely the evolution equations of the metric and of the dilaton field. By 
comparing this ansatz with one made in the previous Section (see Eq. ( p.4|) ) it is evident which 
kind of logic we followed since for r = A = S = Owe exactly recover our Eq. ( p.4| ) . Notice that 
r, A and S are only dependent upon time since we want to accommodate a homogenous dilaton . 
With the ansatz ( |3.8|) we then repeat step by step the procedure outlined in Section 1. First 
of all, we insert Eq. (|3.8|) into Eq. ( p.l4|) and we get the condition which needs to be satisfied in 
order to solve the constraint: 
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By summing and subtracting Eqs. ( p.lO| ) and p.ll| ), we obtain 
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Inserting Eq. (|3.8| ) into Eqs. (|3.1CI| ), we find, respectively: 



A + [A - 0] A + /i tanh/it [2 A - 0] = 0, 

d ( k'\ 3' k' A 
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(3.11) 
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Eqs. (U), (|3lT| ) and (|3A^) can be consistently solved if </)(t) = A(t) = r(t), E = leading to 
the solution 



V\ 
r(t) 



C3 + C2 



jJX 



I ^JL \ C2 C3 /iX 

tanh — H tanh — 
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c| + C2C3 - /ici + Ai V^ T^ '"^ 
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C^ + C2C3 - /iCi + /i^ ux 

2 7 2/i ^^2 



C4 



, ,, (3.13) 

cosh jjLl 

where C3 and C4 are integration constants. Substituting Eqs. (|3.8|) and (|3.13|) into Eq. ( p.9|) and 
Eq. (|2.12|) we obtain, respectively, the following compatibility conditions 



C3 = 3/i^ + C2 — 2/iCi > 0, C4 = C2 — /i^ — 2/iCi > 0. 



(3.14) 



Finally going into Eq. ( |2.13D we deduce that 



2c\ — /i^ + 6/iCi = 3c2 — C3, 



(3.15) 



which turns out to be noting but a trivial linear combination of the two conditions given in Eq. 
( p.l4| ). As expected, Eqs. (|2.10|) and (|2.11|) are then satisfied without imposing any further 
condition on ci, C2, C3 and C4. 

We are now ready to write down the explicit form of our solutions, namely: 



A{x,t) 



^C4gd{pt) 



cosh fit 



sinh ( — 



2 2 

C2+C2C3-A'C1+M 

z^ 



cosh I — 



2 2 



5(a;,t) = e'^^^'^^^*) cosh /it sinh/ix. 



C(x,t) 
0(t) 



-, £2 



cosh /it 

C4 



sinh 



/iX 



e2 + g3 



cosh 



/iX 



c2-g3 



/^ 



5frf(/it) - C5, 



(3.16) 



where gd{fit) = arctan [sinh /it] is the "Gudermannian" |^ or hyperbohc amphtude. Notice that 
the hmit C4 — *> is shghtly ill defined since, from Eq. ( ^.14[ ), C2 might take imaginary values and 



the inequalites are, in this limit, strict equalities. Therefore, the general expression of the solution 
( 3.16| ) also changes, in this limit. 



Since the hyperbolic amplitude always grows, the dilaton will either increase or decrease de- 
pending upon the sign of C4. Therefore, the physical meaning of C5 is connected with the initially 
small (perturbative) value of the coupling constant, whereas C4 is connected with its "kinetic" en- 
ergy. The physical situation we are describing with this solutions is really a regime of very small 
dilaton coupling and in this sense we will fine-tune C5 to a quite small value. In this way higher 
loop (genus) corrections will be automatically subleading. This procedure is very reminiscent of 
what happens in the pre-big-bang scenario |I]] where a (quite long) small coupling regime emerges 
naturally from the tree-level (Kasner-like) solutions. In the pre-big-bang case, the Kasner-like 
solutions have two "branches" separated by a (curvature) singularity, and then, for each growing 
coupling solution there is also a decreasing coupling solution. We notice that, also in this sense, 
our solutions have close analogies to the pre-big-bang case since, depending upon the value of C4 
there are two branches: one of increasing coupling (04 > 0) and the other of decreasing coupling 
(c4 < 0). In the pre-big-bang case the singularity cannot be removed (at tree- level). We will show, 
in the next two sections, that in our solutions the two branches are not analytically connected 
but, at the same time, the curvature invariants, the dilaton kinetic energy are always well defined 
and regular for any x, t at least for a family of solutions contained in Eq. ( ^.ICj ) (see Sections 4 
and 5 for a discussion of this point). 

4 Two physical examples 

From Eq. ( |3.16| ) each particular solution of the system of Eqs. ( |2.4| ) and ( |2.5| ) can be specified 
by fixing ci and C2 and by computing C3 and C4 from Eqs. ( |3.14| ). We tried different possible 
sets of parameters and we found that not all the choices of Ci and C2 will necessarily lead to the 
regularity of the curvature invariants. The class of regular solutions seems to be connected to the 
behaviour of the metric tensor under (discrete) parity transformations (see the following Section). 
At the same time we point out that there are choices of Ci and C2 which lead to completely regular 
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curvature invariants. In this Section we discuss then particular examples of solutions contained 
in Eq. ( |3.1(j| ) that share this quite interesting property. In order to do that we will be inspired by 
some physical requirement, namely we would like to have regularity of the curvature invariants 
and, simultaneously, growing dilaton coupling which corresponds to > 0. 
For the set of parameters given by 



Ci = 6/i, C2 = 5/i, 



C3 



-4/x, C4 = 2v3/i, 



(4.1) 



Eqs. ( p. 14 ) are clearly solved. The explicit form of the solution of Eq. ( 3.16|) can then be written, 
using Eqs. (|4.1|), as: 



C{x,t) 



^2 ^3 gd(^J.t) 



cosh fit 



16 r 



cosh /it 
smh 



cosh 



/iX 



n40 



B{x,t) 



^2 V3 gd{fj.t) 



cosh /it sinh/xx. 



cosh 



fJX 



(p{t) = 2V3 gd{fit) - C5. 



(4.2) 



The dilaton is an increasing function whose maximal value is given by 0(+cx3) = a/Stt — C5 
(recall that \im.t^ -00 9 d{fit) = — (3/2)7r and that limt^+oo gd{fit) = (7r/2)). Now, our discussion 



is based on the (tree-level) action given by Eq. (2J) which is valid provided g{(j)) = exp [0/2] = 
exp [^gd^jjit) — C5] <^ 1. If this is not the case, higher genus corrections should be included. 
Consequently, for the compatibility of the solutions discussed in the present Section with the 
effective description adopted in Section 2 we have to require that the dilaton starts its time 
evolution deep in its perturbative regime and this can be achieved by fixing C5. In the example 
given by the solution (|4.2|) we have that g{+oo) = exp (VSn — C5) <^ 1, which implies that 
C5 ^ \^n. Therefore choosing, for instance C5 ~ IOv^tt we get that g(—oo) = exp (— ISv^vr) and 
(/(+00) = exp (— 9-\/37r). 

Consider now a further example, namely 



Cl 



5/i, C2 = 4/i, C3 = -3/i, 



C4 



VS/i. 



(4.3) 



Using Eqs. ( [4.3|) into Eq. (|3.16|) we get 

l5 



A{x,t) 
C{x,t) 



,V5 gd(tMt) 



cosh fj,t 



cosh nt 

4 r 



cosh ( — 



sinh f — 



n24 



. /iX 

cosh I — 



-\ 7 



B{x,t) = e^ ^'^^^'^ cosh /it sinh/ix, 

0(t) = v^ gd{fit) - C5. (4.4) 
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The calculation of the curvature invariants for the solutions ( [4.2|) and (|4.4| ) can be now per- 
formed. In order to summarise our results we can say that each of the curvature invariants can 
be written (in a generalised notation) as: 

-b 



1 = 12^ /(x,t)[cosh/it]-" 



cosh — 



(4.5) 



(X labels a generic curvature invariant; f{x,t) is a regular expression containing combinations 
of hyperbolic functions and it changes for each specific invariant [see Appendix B]; a and b also 
change depending upon the invariant under consideration, [see also Table 1]). 

From Eq. ( [4.5| ) and from the various forms of f{x, t) which can be deduced from Appendix B, 
we notice that X does not have any pole for any finite value od x and t. Moreover, we point out 
that for large x and t the curvature invariants are suppressed. The magnitude of the suppression 
depends upon a and b which are reported in Table 1 but can also be read-off directly from Eqs. 
( [B.1|) - (|B.4| ). The only nasty thing which could then occur is that for some (large) value of x and/or 
t some of the invariants explode because f{x,t) grows much faster than the suppression factor, 
parameterised, for each invariant, by a and b. Keeping x = Xc (where Xc is some finite value of x) 
and leaving t to run we find that f{xc,t) ~ [cosh/it]^ for large t. Vice- versa keeping t frozen at tc 
(where tc is some finite value of t) we find that (in the worse case) f{x,tc) ~ [cosh^]^. Finally, 
by varying x and t simultaneously we see that f{x,t) can grow, at most like [cosh^]^ [cosh/xt]"^ 
or like [cosh^]^ [cosh/it]^. Therefore we see that the growth of f{x,t) is always suppressed by 
the high (negative) powers (i.e. a and b) of the hyperbolic cosinus appearing in Eq. ( [4.5|) . We 
want now to make few more physical comments concerning the validity of the present solutions. 
The only dimension- full parameter appearing in the solution (|4.2|) and ( [4. 41) is /x = 1/L where 
L is, for the moment, a length scale which completely specifies each solution which describes 
the "spreading" of the curvature invariants around the origin (i.e. x = 0, t = 0). The only 
fundamental scale appearing in Eq. ( p.l|) is the string scale A^. Correspondingly there is also a 
typical energy connected with A^ which is the string energy, namely Es = l/A^. The effective 
action ( |2.1| ) which was the starting point of our considerations holds for energies E <^ Eg. In fact 
for E r^ Eg higher order corrections in the string tension a' = A~^ should be included. In order 
to be consistent with our approach we have to assume that the typical scale L has to be large 
compared to the string scale and this means that L > A^. In the opposite case (i.e. L < Xg), As 
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C4 = 


= 2 v^/i 


RfiuapR^ 


C^.a/3C^'^-'3 


RapR""^ 


i?2 


a 


16 


16 


16 


16 


b 


84 


84 


82 


80 




Ca -- 


= a/5 /i 


-'^fiua/3-'^ 


Cf.uaffC^'''''' 


RalsR'"' 


i?2 


a 


14 


14 


14 


14 


b 


52 


52 


50 


48 



Table 1: We summarize the exponents of the suppression factor appearing in Eq. ( |4.5| ). The 
upper part of the table refers to the case of Eq. ( |4.2| ), whereas the lower part of the table refers 
to the case of Eq. (|4.4| ). The explicit analytical results leading to this table can be found in the 
Appendix B, by setting, respectively, a = c^/ jj, = 2\^ and a = c^j [i^ in Eqs. (|B.1[ - (|B.4|) . 



would not turn out as the minimal length scale. Therefore if L > A^, we will also have X < A~^. 
In Fig. ID the scalar curvature (squared) is reported for the two solutions given in Eqs. ( [4.2|) and 
( [4. 41 ) when L = 10 A^. The corresponding analytical results are reported in Eq. ( [B.4| ) choosing, 
respectively ot = 2^/3 and a = \^. From these analytical expressions we can clearly see that by 
tuning L = l/fi for length scales larger than A^, R^ decreases (in string units) and its spreading 
around the origin increases. We can also look at the behaviour of the time derivative of the dilaton 
which is illustrated ( Fig. ^ and of the related coupling constant (Fig. |^). 

The solutions derived in Section 3 and the examples reported in the present Section can 
be described either in the String frame or in the Einstein frame [^. In the String frame the 
fundamental parameter of the theory is the String length A^. Notice that in this frame the String 
scale is truly a constant, whereas the Planck scale evolves according to a relation which can be 
parameterised (at low energy and low coupling) as Xp{t) = e^Xg. The Einstein frame is defined by 
the conformal transformation which diagonalizes the gravitational action by decoupling the dilaton 
from the Einstein-Hilbert Lagrangian.In the Einstein frame Ap ~ 10~'^'^cm is constant whereas the 
String scale evolves according to As(t) = e^a'Ap (where $ is the Einstein frame dilaton). The 
String frame and Einstein frame dilatons are simply related and, moreover, in four dimensions they 
are indeed equal, whereas the String frame (Qfiu) and Einstein frame (G^jy) metrics are connected 
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Figure 1: We plot the R^ invariants corresponding to the solution (|4.2| ) (left) and to the solution 
( [4.4|) (right). The analytic form of these invariants is given in Appendix B (Eq. ( p.4|) ) by choosing 
a = 2-\/3 (for the plot at the left) and a = a/5 (for the plot at the right). Here we also took 
L = 10 A, 
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Figure 2: We plot the dilaton kinetic term for the two exact solutions discussed in this Section. 
The left figure corresponds to the case given in Eq. ( [4. 21 ), whereas the right one corresponds to 
Eq. ( ^4.4| ) . We notice the absence of any singularity. 
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Figure 3: We report the evolution of the dilaton couphng in the solutions given by Eq. (|4.2| ) (left) 
and by Eq. (|^ (right). We took C5 = 10 ^3 tt (left) and C5 = 10 ^5 tt (right). Provided the 
dilaton coupling starts its evolution deep in the perturbative regime, g{(j)) ^ 1 for every time. 



by a conformal rescaling involving the dilaton coupling {g^i, = e'^G 



' l^uj 



We transformed our solutions in the Einstein frame and we investigated the emerging picture 
which can be summarized by saying that the curvature invariants are also regular and everywhere 
defined. At this point, we would be tempted to call the solutions we just discussed "singularity- 
free" . In order to name a generic solution of the low energy String effective action "singularity- 
free" we should discuss the geodesic completeness of the metric both in the Einstein and in the 
String frame pictures. From the point of view of the interpretation of our solutions, moreover, 
the geodesic completeness in the String frame seems even more crucial, since strings moving in 
curved space follow geodesies in the String frame but not in the Einstein frame ||2^. To show 
that the geodesies do not hit a boundary in both frames is anyway a quite long computational 
task. We prefer to postpone it to further investigations mainly because the calculations reported 
here are already quite heavy. We do not have, at the moment, any argument either favouring or 
disfavouring the geodesic completeness and we believe that this can only come out from a precise 
calculation. 

There is a second very important issue. The Einstein frame picture of our solutions involves 
an energy momentum tensor which does not violate any of the energy conditions which are part of 
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the hypotheses of the powerful Hawking- Penrose (HP) theorems P5| , |26| . Therefore the question 
which naturally arises is which other hypothesis of the HP theorems is violated. There are indeed 
examples of truly non-singular (inhomogeneous) solutions in general relativity which have an 
energy-momentum tensor of a perfect fluid I^Tj. In that context it was rigourously proven that 
the so-called initial or boundary condition of the HP theorems was the one failing to hold. This 
would mean, in our context, that we have to fine-tune (as we had!) the initial conditions quite 
heavily in order not to bump into a singularity. This point is indeed closely connected with the 
behaviour of the geodesies and with their completeness and therefore it certainly requires further 
studies. The same considerations should be performed in the String frame which is perhaps the 
most useful for the physical intuition (as it was noticed in the case of the pre-big-bang scenario |jl|). 
We are not aware, at the moment, of any general singularity theorem in this case and therefore 
we are not in condition of making any kind of general statement going beyond the solutions we 
just presented. 

5 Parity invariant solutions with growing coupling 

The two physical examples presented in the previous Section share a quite amusing physical 
property, namely the invariance of the line element corresponding to each one of the solutions 
given in Eqs. ([4.2|) and ( [4.4|) under discrete parity transformations. In fact by computing, from 
Eq. ( p.7|) , the line element associated to the examples reported in Eq. ( [4.2|) and (|4.4|) we find 



that it is invariant under the transformation x —>■ —x. Surprisingly enough, we will show that 
parity invariant metrics lead to regular curvature invariants. This will allow us to define a class 
of regular, inhomogeneous models with growing coupling. 



We start by noting that, from the structure of the metric written in Eqs. (|3.16| ) the diagonal 
elements of the metric tensor given in Eq. ( |2.6| ) are even under parity transformations provided 
the spatial (x) dependence appears in functions which are themselves even under parity. Given 



the form of the metric (p.l6|) we see that, in order to preserve parity, the spatial dependence 



has to occur either through powers of the hyperbolic cosinus (which is parity even) or through 
even powers of the hyperbolic sinus. We then impose the following two extra conditions to the 
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parameters Ci, C2, ,03, C4 specifying the solutions (|3.16|) : 

C2 + C2C3 - /iCi + /i^ = 0, C2 + C3 = /X (5.1) 

In Eqs. ( J.16D the first condition of Eq. ( |5.1| ) sets to zero the power of the hyperbohc sinus 
appearing in A{x, t) whereas the second condition sets to one the power of the hyperbohc cosinus 
appearing in C(x, t). Notice that with this choices the components of the metric are all even under 
X -^ —X. Therefore, the full system of conditions obeyed by the solution given in Eq. (|3.16|) is 

(5.2) 



^^^'=3+r^V-2^^^ 



C4\2 



Cl 

/i 



/^ 



C2WC3 



(5.3) 
(5.4) 

, (5.5) 

Now calling, for simplicity a = 04//^, 7 = C2//i, S = Ci//i and and ehminating C3//X through Eq. 
(^.2|) from the previous equations we obtain the following system: 



^1 + r-i) = 1. 



7 + 1 
,2 



a 

27 + 2 






26, 



26. 



(5.6) 
(5.7) 
(5.8) 



Eqs. ( |5.6| ) and ( |5.8| ) are clearly not independent and, therefore, the solution given in Eq. ( |3.16| ) can 
be expressed only in terms of one (physical) parameter which we choose to be a. This turns out 
to be a physically interesting parametrization since /i^a^ = 0^(0). In this sense |a| estimates the 
(maximal) energy density of the dilaton background in string units. The sign of a automatically 
selects solutions either with growing coupling or with decreasing coupling. Thus, using Eqs. 
( p. 61) and ( |5.7| ), we have that the metric of Eq. ( p.l6| ) can be re- written as 



A{x,t) 
B{x,t) 
C{x,t) 



3±« gd(Ati) 



cosh fit 



2±Vc?+4 r 



cosh I — 



-1 2Va^+4[Va^+4±l] 



e±"srf(A«t) cosh fit sinh/ix, 

cosh fit sinh 

±a gd{fit) — (3. 



fix 



cosh 



fix 



l±2v'a2+4 



(5.9) 
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(the plus and minus signs correspond, respectively, to solutions with growing and decreasing 
dilaton coupling). In Appendix B we computed the curvature invariants in the case of growing 
coupling solutions. We indeed found that the requirement of parity invariance of the solutions 
seem to be winning since the curvature invariants are regular. Notice that the curvature invariants 
of the two previous examples (reported in Section 4) can be derived from the general result 
of this Section by setting, respectively, ot = 2\^ and a = ^/E. Also in the case of the class 
of solutions discussed in the present Section the curvature invariants vanish asymptotically as 
can be argued from the results reported in Appendix B. The curvature invariants can be also 
computed in the case of decreasing dilaton solution. We did this exercice and we found, again 
regular solutions. From the more mathematical results reported in this Section (and from the two 
examples of Section 4 illustrating the physical properties of our solutions) we see that something 
quite unexpected happens, namely the fact that the regularity of the curvature invariants can 
already be achieved at tree level (without any curvature correction). Moreover, as shown in 
Section 4 for two particular cases, we can set initial conditions in such a way that the coupling 
constant is always much smaller than one. We notice that these results are not in contrast with 
previous discussions |^ suggesting that the only way of regularizing the curvature invariants was 
to invoke higher curvature corrections. In fact those studies were always postulating a completely 
homogeneous (and often isotropic) four dimensional geometry. The physical picture emerging from 
our detailed analytical calculations is that the inhomogeneities are quite an essential ingredient 
in order to regularize the curvature invariants. On a purely physical ground, moreover, it does 
not seem unreasonable to have strong inhomogeneities at very small scales. As we noted in the 
introduction our results confirm and extend previous ideas concerning possible avoidance of the 
singularities in completely inhomogenous string cosmological models ||12||. In agreement with |12 



we showed that regular solutions are possible at tree-level. At the same time we showed that this 
conclusion holds also in the presence of a growing coupling solution and in the absence of any 
extra field. We believe that our new results should be taken as a serious motivation in order to 
perform a fully consistent numerical study of inhomogeneous string cosmological models. These 
studies might extend (and perhaps change) some of the conclusions drawn in the completely 
homogeneous case [^]. In particular, numerical studies might clarify if it is at all possible to 
find inhomogeneous examples evolving at late times towards a completely homogeneous state. 



18 



In fact a quite unphysical feature shared by our models and by the ones of |jT2[ is their eternal 



inhomogeneity. It is certainly tue |]T^ that the homogeneous limit of a inhomogeneous model is 
already quite delicate in general relativity. In spite of the possible subtleties we can say that a 
good measure of how inhomogeneous and anisotropic is a solution at late times can be the ratio 
of the Weyl invariant {C^^^^C^uap) over the Riemann invariant {R^'^"'^ R^i^ap) ■ In our case, from 
the results of Appendix B, this quantity goes to a constant at for t —>■ +00. A similar behaviour 
can be found in []T2[. Numerical studies could also give us some clue concerning the role played 
by the compactification radii which we took frozen but which should be included in the game of 
inhomogeneous models. Up to now this possibility was never considered. 

6 Concluding remarks 

In this paper we gave some examples of regular dilaton driven solutions of the low energy beta func- 
tions without including any dilaton potential, any higher a' (curvature) corrections, any higher 
genus corrections, any antisymmetric tensor field any dilaton potential. We took the six com- 
pactification radii to be constant and we looked at some special class of inhomogeneous metric. 
We showed that if the metric is sufficiently inhomogeneous the curvature invariants are bounded 
for every x and t. The only definite conclusion that we feel like stating is that regular (but in- 
homogeneous) solutions of the evolution equations of the dilaton and of the geometry are more 
common (at tree-level) than in the completely homogenous case where no-go theorems probably 
apply. Since in our examples the coupling is always very small and the curvature invariants always 
smaller than one (in string units) the tree-level treatment seems, a posteriori, not meaningless. 
The curvature invariants (including the Weyl invariants) vanish, asymptotically for large times. 
Moreover the metric is invariant under parity transformations. 

We leave to future studies many unanswered questions. Can we promote the regular solutions 
we spotted to true singularity free solutions? If this is the case are the geodesies complete in the 
Einstein frame? Which hypotheses of the singularity theorems turn out to be violated? Can we 
formulate the singularity theorems (in a consistent way) directly in the String frame or we have 
always to go to the Einstein frame? Should we consider regular only those solutions which are 
bounded in both frames? 
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APPENDIX 

A Christoffel symbols and Ricci tensors 

In this Appendix we will give the general formulas for the Christoffel symbols and Ricci tensors 
computed in the case of the diagonal metric given in Eq. (|2.6|) with two commuting space-like 
killing vectors which are hypersurface orthogonal. 

The Christoffel symbols can be easily derived from the metric ( ^.61 ) and they are reported in 
Table |^, whereas the Ricci tensors and the curvature scalar can be obtained, after some algebra, 
once the Christoffel symbols are known. Our conventions for the Ricci tensor will be -R^i, = Raav°' 
(where R^^a ^ = '^ti^ua^ ~ ■■■ is the Riemann tensor). The (xx), {yy), (zz) and (00) components 
of the Ricci tensors are, respectively 

1 f .^N2 ^B A ^A\^ A'B' ^B\^ ,C\^ A" B"\ ,^ ^, 

y~~YA\BC~^c' ^B^c~^~c^^~c' ~^~~cy ^ ' 

1 fC^_BC B _C B'C'_C' B" C"\ 
^^~ 2AY& BC^B C^BC ^C' 5 + C j' ^ '^^ 

«S = J^i'l' + AB + 'b> - 'c' -a-'b- (t' + A ¥ + -a]' <*-^' 

The (Ox) component and the curvature scalar will be instead: 

o_J_iBA AB^ B^B_CC_ B^) 
^^~2A\BA^ AB^ BB CC 5 J' ^^'^^ 

^^i^^^^{2(^)^ + (|)-(§)-2f4|-4)^-(f)^.(^)^.2f.4^ 
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■px _ 1 9 log A 
'- XX 2 dx 


yx _ C Bdlog{B C) 

yy 2 A ax 


J.X _ B aiog(B C) 
zz 2 A C ax 


-px 1 9 log A 
^ xQ 2 a t 


-px 19 log A 
^00-2 ax 


p^ _ 1 log (B C) 
^ xy 2 ax 


py _ 1 log (B C) 

'- yO 2 at 


-pz _ ld\og{B/C) 
^ xz 2 ax 


J.Z _ laiogiB/C) 

^ zo 2 at 


pO _ 1 <9 log A 
XX 2 dt 


pO _ C B51og{C B) 
^ yy 2 A at 


pO B a\og{B/C) 
'- zz 2 A C dt 


pO 19 log A 
'- xO 2 dx 


pO 19 log A 
■^00 2 9t 



Table 2: We report the Christoffel symbols computed from the metric given in Eq. (|2.6| ). 

B Curvature invariants for growing dilaton solutions 



From the (parity invariant) class of solutions given in Eq. (|5.9| ) we can compute the curvature 
invariants, namely curvature scalar together with the squares of the Riemann, Weyl and Ricci 
tensors. The result is the following 



R.uapR^""^ = ^e-^"^'^^''*) [cosh /it] 



-8-2Va^+4 



cosh 



fix 



-4-4Va2+4(l+Va2+i) 



si[a 






RaffR"' 



+ S2{a) cosh 2 fit + S3 (a) cosh4/it + ^4(0;) cosh/ix + 35(0;) cosh2/ia; 

+ Sei^a) cosh 2 fit cosh 2 fix + SY{a) cosh 2 fit cosh fix 

+ Ssia) cosh 4:fit cosh fix + sg{a) sinh fit + sio{a) sinh fit cosh 2 fix 

+ sii (a) sinh /it cosh /ix > , 



t!_ ^-2agd(pt) 

96 



e -"^"^f-"' [cosh /it] ^ 



-2Va2+4 



cosh 



/iX 



-4-4\/c?+4(l+\/c?+4) 



wilal 



+ 1^2(0) cosh4/it + ii'4(a) cosh/ix + w^i^a) cosh2/ix 
+ Welct) cosh 2 fit cosh 2/ix + ^7(0) cosh 2 fit cosh /ix 
+ Ws{a) cosh 4/it cosh /ix > , 

-2-4Va2+4{l+V«2+4) 



„2,,4 
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[cosh fit] 



-8-2Va^+4 



cosh — 



(B.l) 



(B.2) 



ri(a) 



+ r2(a;) cosh2/it + r^i^a) cosh/ix + r4(a) cosh 2/it cosh /ix 
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R' 



+ r^{a) sinh fit + r^i^a) sinh fit cosh fix 



4 4 

a jjL e 



cosh 



fix 



-4Vc?+4(l+v'c?+4) 



(B.3) 
(B.4) 



The a-dependent coefficients appearing in our expressions are reported in Table |^. Notice that 
the two examples specifically discussed in Section 4 arise by setting, respectively, a = 2-\/3 and 
a = ^/5. Plugging these values in Eqs. (p.l| )- (|B.4|) we get the corresponding curvature invariants. 
We want to mention that we did the calculations of the curvature invariants also in the Einstein 
frame after having transformed the solutions according to the general conventions (discussed in 
Section 4) which relate the two frames. The curvature invariants turned out to be regular also in 
the Einstein frame. 
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si{a) 


1152 + 416a^ + 824a^ + Ga"^ + V«2 ^ 4(594 + gga^ + Ua^) 


S2{a) 


384 + 536a' + 124a^ + Sa*^ + Va' + 4(96 + 164a2 + 20a^) 


S3{a) 


384 + 216a2 + 38a^ + 2a'' + V^^ + 4(168 + 66^2 + 6a*) 


s^ia) 


48 - 108a^ - 78a* - 6a'* - Va^ + 4(24 + igOa^ + 18a*) 


S5(«) 


192 + 104a2 + 16a* + y/a' + 4(192 + 36a2 + 4a*) 


seia) 


72a2 + 12a* + Va' + 4:{36a' + 4a*) 


S7{a) 


-1728 - 688a2 - 120a* - 8a^ - (864 + 232a2 + 24a*) V«2 + 4) 


ss{a) 


-144 - 108a2 - 18a* - (72 + 42a2 + 6a*) Va^ + 4 


sg{a) 


96a^ + 24aVtt2 + 4 


sio(a) 


32a^ + 8aVa2 + 4 


sii(a) 


128a^ + 32a V«2 + 4 


Wi{a) 


3456 + 1488a2 + 216a* + 12a'^ + Va' + 4(1512 + 390a2 + 30a*) 


W2{a) 


1152 + 1128a2 + 282a* + 18a^ + Va^ + 4(288 + 300a2 + 48a*) 


W3{a) 


1152 + 648a2 + 114a* + 6a'' + y/a' + 4(504 + 198a2 + 18a*) 


W4{a) 


-144 - 492a2 - 266a* - 180^^ - (72 + 378a2 + 54a*) Va^ + 4 


W5{a) 


1152 + 552a2 + 94a* + 6a'* + Va^ + 4(576 + 204a* + 24a*) 


weia) 


216a^ + 74a* + 6a'* + \/a^ + 4(108a2 + 24a*) 


W7{a) 


-5184 - 1824a2 - 1128a* - 24a'* - Va^ + 4(2592 + 888a2 + 72a*) 


ws{a) 


-(216 + 126a2 + 18a*) Va2 + 4 


ri{a) 


22a2 + 2a* + 4a Va^ + 4 


r2{a) 


80 + 22a2 + 2a* + (32 + 4a2) Va^ + 4 


rsia) 


-80 - 14a^ - 2a* - (32 + 4a2) Va^ + 4 


Tiia) 


-14a2-2a*-4aV«2 + 4 


r5{a) 


32a + 8aVa2 + 4 


raia) 


32a + 8aVa2 + 4 



Table 3: We report the coefficients appearing in the expression of the Riemann, Weyl and Ricci 
invariants given, respectively, in Eqs. (p.l| ), ( p.2|) and ( p.3|) . 
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